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Abstract 

Chaos in Robertson- Walker cosmological models where gravity is coupled to one 
or more scalar fields has been studied by a few authors, mostly using numerical 
simulations. In this paper we begin a systematic study of the analytical aspect. 
We consider one conformally coupled scalar field and, using the fact that the 
model is integrable when the field is massless, we show in detail how homoclinic 
chaos arises for nonzero masses using a perturbative method. 
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I. Introduction 

In Ref. 1, Calzetta and El Hasi presented an argument, supported by numerical evidence, 
for the appearance of chaotic behaviour in a spatially closed Robertson- Walker (RW) 
cosmological model, filled with a conformally coupled massive scalar field. Because cosmo- 
logical models of this type are practically the simplest possible ones, but give in some sense 
a rough, smoothed out indication of the average properties of a more realistic cosmology, 
they have frequently been used as a first testing ground on which to try out approaches 
to both classical and quantum problems in cosmology. In addition, because of their kine- 
matical simplicity, one would expect the dynamics of these models to be equally simple. It 
thus may come somewhat as a surprise that chaos should appear at this level, although it 
is known that it can appear in a nonlinear dynamical system with two degrees of fredom, 
and one may suspect that it is in fact generic for relativistic systems. 

The study of chaos in relativistic cosmology began more than twenty years ago with the 
study by the Russian school 2 and Misner's group 3 on diagonal Bianchi type IX (mixmaster) 
vacuum cosmologies. This research is, by itself, very interesting and still active, 4 but the 
hopes initially placed in these models, essentially the solution of the horizon problem 
from the mixing effect associated with chaos, could not be sustained, and the Bianchi 
IX cosmologies have been mostly used as toy models, rather than first approximations 
from which to extract estimates of physical quantities. On the other hand, although the 
question of which RW model best approximates our inhomogeneous, anisotropic universe 
is not yet settled, there is agreement on the relevance of RW models both for classical and 
quantum cosmology, 5 and physical predictions are often made based on them. Therefore 
the presence of chaos, in this kind of cosmologies, enlarges the avenues of research in 
fundamental subjects such as the arrow of time and the issue of irreversibility in the 
evolution of the universe, or the transition from the quantum regime to the classical one, 
with the corresponding appearance of decoherence and correlations and related phenomena 
like particle creation, and has been used in considerations related to physical phenomena 
in inflationary universes. 6 ' 7 ' 8 

In their work, Calzetta and El Hasi used mostly a numerical method to investigate 
the chaotic behaviour of the model, but also sketched an analytical, perturbative method 
by which they considered the model as a perturbation of the integrable one obtained with 
a vanishing mass for the scalar field, and split the Hamiltonian of the system into an 
integrable part and a coupling term. With a perturbative argument commonly used in the 
treatment of near-integrable systems, based on the Melnikov criterion for homoclinic chaos 
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and Chirikov's resonance overlap criterion, 9 they argued qualitatively that the KAM tori 
of the integrable part are destroyed and replaced by stochastic layers under the effect of 
the perturbation. The resulting strong indications of chaotic behaviour were confirmed by 
the numerical analysis, which allowed them to go beyond perturbation theory. 

The generic transition to chaos in a perturbed integrable Hamiltonian system ( "soft 
chaos") can be seen as taking place in two steps. 9 ' 10 In the first one, a perturbation of 
the Hamiltonian, with a nonvanishing Fourier component in resonance with one of the 
tori of the unperturbed system, modifies the dynamics in a neighborhood of this torus to 
give a new integrable system, which however has heteroclinic orbits even if before there 
were none. These orbits are separatrices, which act as seeds of chaos when in the second 
step one takes into account the effect of the remaining components of the perturbation 
on the new integrable dynamics; it is here that one may use the Melnikov method to 
determine whether in effect the dynamics becomes chaotic. What this method provides is 
a topological criterion for detecting a chaotic "tangle" in the orbits near the separatrices, 
and is therefore free of the coordinate ambiguities that have arisen, e.g., in the study of 
the chaotic nature of Bianchi models by other methods; for a quick introduction to the 
Melnikov method see, e.g., Refs. 11 or 12, and Ref. 13 for a more technical exposition. 

The goal of this paper is to continue the work in Ref. 1 by starting to examine in 
detail the analytical treatment of the model; in future work, we will systematically extend 
our work to more general models. We will use a different perturbative expansion of the 
Hamiltonian, and show the onset of chaotic behavior by an explicit computation of the 
Melnikov integral. Other recent work on chaos in Robertson- Walker models coupled to 
scalar fields has focused on their application to inflation, again using mostly numerical 
methods, and an analysis based on the Painleve analysis of differential equations, 14 which 
confirmed the non-integability of the models. The main advantage of our work is the 
amount of information it can give us regarding the resonances responsible for the onset of 
chaos and the chaotic region of phase space. 

The paper is organized as follows: in section II, we set up the perturbative method 
we will use, by choosing appropriate variables to describe the cosmological model and sep- 
arating the corresponding Hamiltonian into an "unperturbed" term and a "perturbation;" 
in section III we analyze the first step above for our example, i.e., we discuss the local 
dynamics in the vicinity of one of the resonant tori, while in section IV we develop the 
second step, discussing the effect of other perturbation terms on the local dynamics, which 
becomes chaotic as Melnikov's method shows; section V contains concluding remarks on 
possible future developments. 
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II. Dynamical variables and Hamiltonian: 
setup for the perturbative treatment 

Let us consider the closed Robertson- Walker (RW) metric for a homogeneous and isotropic 
universe, written in the form 

ds 2 = a 2 {t) [-dt 2 + d X 2 + sin 2 X (d6> 2 + sin 2 6> dip 2 )} , (2.1) 

where 0<</?<27r,0<#<7r,0<x<7r are the angular coordinates on S 3 , and t is 
the usual "conformal time" (see, e.g., Ref. 15). We shall consider only models possessing 
a cosmic singularity a = 0, chosen for convenience to occur at t = 0. We shall also assume 
that after the big crunch a new cosmological cycle begins, but a changes sign, in such a 
way that the evolution is smooth, and this process is repeated an infinite number of times. 
This well-known possibility of extending the evolution beyond the a = singularity is an 
artifact of the symmetry possessed by the RW models (we have no grounds for believing 
that it can be generalized), but it is convenient for us because it allows us to assume that 
time runs up to t = +oo, and thus to use techniques for detecting chaos, including its very 
definition, that would otherwise be inapplicable. Any conclusion we draw from the model 
would of course not have to rely on this extension; the time scale for any chaotic behavior 
to set in must be shorter than one cycle for its physical consequences to be meaningful. 
We start with the Einstein-Hilbert action for the gravitational field, 

and the action 

S^g\ = -\ Jd*x y /=g[g ab V a *V b *+(/JL 2 + iR)&] (2.3) 

for a real, conformally coupled scalar field of mass fx. To write down the Hamiltonian, use 
the fact that the square root of the determinant of the metric (2.1) and its scalar curvature 
are, respectively, 

y/=g = a 4 sin 2 x sin6> , R = 6 (j^ + , (2-4) 

where an overdot denotes a t-derivative, and reparametrize the scalar field (which must 
be homogeneous for consistency) by $ i— > </> := y^4nG/3a^. Then we obtain, up to a 
constant overall factor, 

H(a, 0; 7T, p) = \ [ - (tt 2 + a 2 ) + (p 2 + <j> 2 ) + ^a 2 ^ 2 ] , (2.5) 



4 



where n and p are the momenta conjugate to a and 0, respectively. Here, as usual with 
general relativistic systems, time reparametrization invariance requires that we impose 
that all solutions satisfy the constraint 

H = , (2.6) 

which is just the (0, 0)-component of the Einstein equation. 

The system described by the Hamiltonian (2.5) consists of two harmonic oscillators 
(of which one is "inverted," as is to be expected from any degree of freedom related to 
the spatial volume element), coupled through a term proportional to p 2 . For p 2 = 0, 
the system is trivially integrable; our main goal is to analyze the effect of the coupling 
term, at least for small values of p 2 . However, in view of the perturbative treatment to 
be carried out, the form (2.5) has two disadvantages: the "perturbation" term p 2 a 2 (p 2 is 
not always small, in the sense that for a fixed value of fx 2 it cannot be uniformly bounded 
for all orbits of the unperturbed system; and the latter system is degenerate, in the sense 
that the frequencies of both oscillators are constant. We will not attempt to deal directly 
with the first inconvenience, which at any rate is not really a problem for perturbations 
of a fixed orbit; we will instead perform a change of variables that will allow us to isolate 
a different integrable Hamiltonian from H, with enough of the degeneracy removed for 
us to be able to use the Melnikov method of detecting homoclinic chaos. Regarding the 
perturbation, it will be sufficient for us to know that it is a well-behaved function on each 
trajectory. 

Following Ref. 1, we first replace the dynamical variables p and <fi by new variables j 
and <p, respectively, defined as 

4> = \ / — sin ip , p = a/2 uj cos (p , (2-7) 
V ui 

where u = a/1 + p 2 a 2 is the instantaneous frequency of the field; thus, this is a first 
approximation to action-angle variables for the scalar field, which takes into account the 
coupling term, but not the fact that a is a dynamical variable, and it is a key trick in this 
treatment of the model. To complete the canonical transformation, we must introduce a 
new momentum variable conjugate to a, given by 

2 

P:=7r-|^sin2^, (2.8) 

and arising from the fact that u is a-dependent. Now the Hamiltonian in terms of the new 
variables can be written as a sum H = H Q + 5H, of an unperturbed Hamiltonian 

H (a,cp;P,j) = -\ (P 2 + a 2 ) +j^/l + ^ a 2 , (2.9) 
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which is obviously integrable, since it has the two commuting constants of the motion H Q 
and j, and a perturbation 

2 

(l-cos4^). (2.10) 

However, the system (2.9) is not easy to integrate; since we will need to perform explicit 
calculations, and in particular to write it in terms of action-angle variables, we find it 
convenient to simplify it somewhat. 

In Ref. 1, this simplification was achieved by noticing that for most of the unperturbed 
motion, the value of u is very large, and using the u ^> 1 approximation to H Q . In addition, 
because the treatment was a perturbative one in the parameter ,u 2 , or more precisely fj, 2 j, 
the second term in the Hamiltonian perturbation (2.10) was dropped, since it contains 
higher powers of /j, 2 than the first one. Notice that, implicit in this way of proceeding is the 
idea that one considers the in H Q as a fixed parameter, while in 5H, y 2 =: e is rendered 
temporarily independent of the first one for the sake of the perturbative treatment, and 
allowed to vary in a neighborhood of fx 2 = e = 0. (One may be able to set up a perturbation 
theory for Hamiltonians of the type H Q {q,p;e) + e 5H{q,p), where H (q 7 p;e) is integrable 
for all £, but to our knowledge such a theory is not available yet.) 

Here, we will follow a different path. Since we are going to expand the perturbation in 
powers of \i 2 and neglect terms higher than fi 4 , we will do the same for H Q , and retain there 
only the terms up to order /i 2 while we will put in the perturbation the higher order ones. 
The new unperturbed system obtained in this way will be explicitly integrable. As the 
perturbation parameter, we will use /U 2 , rather than the dimensionless n 2 j; this is mainly to 
avoid having to explain already at this stage what values of the variable j we are interested 
in. At any rate, strictly speaking we will not have to talk about s — /j 2 being "small" (if we 
did, we might say that /x/%> <C 1, where fi P is the Planck mass); the perturbative results 
we will obtain are valid for e in some interval near zero, and such topological statements 
are "dimensionless." Expanding a/T+ \x 2 a 2 = 1 + \\j?a 2 — \\i 4 a 4 + o(p 4 ), we obtain a 
new breakup of the Hamiltonian into 

H = H Q + 5H, (2.11) 

where the unperturbed part is now 

H (a, <p; P, j) = -\ (P 2 + a 2 ) + j (1 + \ fa 2 ) , (2.12) 

and the perturbation differs from (2.10) by terms proportional to /U 4 a 4 and of higher order 
in /I 2 , namely 

5H = -\n 2 aPj sin2v?+/ [\ a 3 Pj sin 2^ - \ a 4 j - ^ a 2 j 2 (1 - cos4ip)]+o(/j 4 ) . (2.13) 
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5H(a 1 <p;P l3 ) 



fx z aPj 
'2(l + fi 2 a 2 ) 



sin 2cp — 



4(l + fi 2 a 2 ] 



One of the action variables of the new unperturbed system is j itself; to find the 
other one, we compute the integral k := (2tt)~ 1 f P da, along the trajectory specified by 
the values of the conserved quantities j and H = const =: h (ultimately, we will be 
interested only in trajectories near h = 0, in order to satisfy the constraint H + 8H = 0, 
but for now we keep the treatment more general), for which 

P(a) = ±^2 (j - ho) - (1 - » 2 j) a 2 . (2.14) 

The integration yields 

k = ~j=^= , (2-15) 
and the angle variables canonically conjugate to k and j, respectively, are given by 

9 = arctan (V 1 ~ n 2 j a/P) (2.16) 

as can be checked by direct computation of the Poisson brackets. Notice that \i 2 here is 
not the infinitesimal e, but we do want to consider it as being small, so for j not too big 
1 — fi 2 j > 0. With these new variables we can now write down our final expression for the 
unperturbed Hamiltonian, 

H (k,j)=j-ky/l-l* 2 3 • (2-18) 
Choosing a value for h Q will then impose a relationship between j and k; in particular, for 
h 0, since j is positive — see (2.12) — k must be positive as well. For the perturbation, 
we invert the transformation (2.15)-(2.17), and substitute in (2.13), which, keeping only 
the lowest order terms in e = [i 2 , gives 

5H = \ e kj [cos(2# + 26) - cos(2# - 25)] + \ e 2 [-§ k 2 j - \ kj 2 + 

+ (2 k 2 j + \ kj 2 ) cos 26 - \ k 2 j cos 49 - \ k 2 j cos(2o - tt/2) + \ kj 2 cos4o + 

+ k 2 j cos(2^ - 28) - k 2 j cos(2# + 28) - | kj 2 cos(2^ + 48) - | kj 2 cos(2# - 48) + 

+ *f k 2 j cos(40 + 28 + if;) - ^ k 2 j cos(4# -28 + ip)] + 0(e 3 ) , (2.19) 

where ip := — arcsin(l/A/5), and we have chosen to write the various terms in the form 
that will be most useful in the following. This completes the setup for the perturbative 
treatment of the model. The dynamics of the unperturbed Hamiltonian H is trivial, and 
gives a conditionally periodic motion 9 = 9 + u k t, 8 = 8 + u-t, with frequencies given by 



d_Ho 

dk 

dH Q = v 2 k + 2^/1-^j 
dj " 2 v 'l //-'./ 
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(2.20) 



In the next section we will start to take a look at the perturbed dynamics. 



III. Effect of the resonant terms in the perturbation 

The effect of the perturbation SH is felt in particular on the rational tori of the unperturbed 
dynamics, the ones where the motion becomes periodic, because the frequencies (2.20) 
associated with the two degrees of freedom are related by the resonance condition 

for some integer numbers n and m Q . Using the explicit form of the frequencies in the 
resonance condition, we obtain a relationship between the action variables k and j, namely 
that 

2(l-/i 2 j) _ m 



fi 2 k + 2 v 7 ! - ^3 n o 
must be a rational number. The latter equation can be expressed as 



(3.2) 



= 2 IhL (i _ ^j) _ 2 v/T^j , (3.3) 



m 



and, since k must be positive in order for h 0, we see that a necessary condition is that 
n and m have the same sign — we can take them to be positive — , with n strictly greater 
than m . 

To find out which of the resonant tori of the unperturbed system are affected by 5H, 
we must consider the perturbation (2.19) as a Fourier series with respect to the angular 
variables 9 and 5, of the form 

oo 

8H(6,8;k,j) = J2tH^(e,8;k,j) 

6H®(9, 6; kj) = e*J2 V<&(k,j) cos(n# + mS + , (3.4) 

n,m 

and check which of the Fourier coefficients Vnm(k,j) are non-vanishing. The order e terms 
in (2.19) are 

5H (1) = \ekj [cos(2# + 26) - cos(2# - 25)1 . (3.5) 

Therefore, to order [i 2 in the perturbation, the resonant tori which are broken by the 
perturbation are those for which 

n = ± m = 2 . 

However, neither of these pairs satisfies the conditions given below (3.3), which means that 
these resonant tori are not present in the physically relevant region near h = 0. 
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This forces us to consider the Fourier components of order e 2 in the perturbation 
(2.19). In order to look at their effect on the dynamics, however, we must somehow consider 
them as perturbations of a system which includes the order e terms. The dynamical 
behavior of such a system is qualitatively similar to that of the unperturbed one in the 
region near h = 0, because H has no tori there which resonate with the order e terms, 
and to set up the system we can use the normal techniques for non-resonant perturbations 
of integrable systems. 

The method (see, e.g., Ref. 16, §5.10) consists in performing an e-dependent canonical 
transformation to new coordinates (6',5';k',j'), in which the new Hamiltonian has no 
perturbation terms of order e, and is obtained by using k'6+j'5+&(6, 5; k',j') as generating 
function, where 

*(*. * - - E ^f^g^T-Tj Mr* + mi + . (3.6) 

While the new canonical coordinates cannot be explicitly expressed in terms of the old 
ones, they differ by a series of powers of e, and for our purposes it will be sufficient to 
calculate the term of order e. A calculation shows that the new unperturbed Hamiltonian 
H' G (k',j') = H Q (k',j'), i.e., it is given by the same function (2.18) as in the previous 
variables, while the perturbation is now of the form 



5H' 



£ 



V^\k\f) + V£\k',j') cos 26' + V£\k',j') cos(25' - tt/2) + 

+ V^\k',j') cos 49' + V^\k'j') cos 45' + 

+ V22 2 \k'j') cos(26>' + 25') + v£l{k',j') cos(26>' - 25') + 

+ V 2 f\k'j') 008(2(9' + 45') + V^l(k'J') cos(2#' - 45') + 

+ V^\k',j') cos(40' + 25' + VO + V^ik'j') cos(4#' - 25' + </>) + 

+ V^\k',f) cos(40' + 45') + V^l(k',f) cos(40' - 45')} + 0(e 3 ) . (3.7) 



Since the unperturbed Hamiltonian is the same as before, the condition for a resonance at 
h ~ is still n > m > 0, and the only term in (3.7) satisfying this condition has 

n = 4 , m = 2 , 

in which case 



tl 2 k Q = 2 y/1 - ^J Q (2 y/1 - - 1) , (3.8) 
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is positive if fj> 2 j < 3/4. We thus choose to study the motion near the resonant torus 
with action variables (j ,k Q ) characterized by the integer numbers (n ,m ) = (4,2), and 
for which the perturbation coefficient is 

v£\k',j') = vS\k',j') = 4k' 2 j'. (3.9) 

(It is not necessary for us to write down explicitly all the functions Vmn (k',j') here, but 
some of them coincide with the Vmn(k',j') in (2.19).) As far as this local dynamics is 
concerned, then, we could have used directly the Hamiltonian H (k,j) with perturbation 
, j) as in (2.19), without worrying about the presence of the e term. From now 
on, we will drop the primes on the variables just introduced. 

Since the perturbation with coefficient (3.9) is resonant, it does change the local 
dynamics qualitatively near the chosen torus. Following the standard procedure for such a 
case (see, e.g. Ref. 9), we go over to a set of canonical variables (7, 5;K, J) adapted to the 
resonant torus, chosen so that one of the momenta will still be a constant of the motion 
under the resonant term in the perturbation: 
k — k 

K := ^ = \ (k - k ) 7 := n 6 + m 6 + % = 4 9 + 2 5 + % 

m (3 - 10) 
n 

where t(j Q is some arbitrary, fixed angle. 

If we suppose that K <C 1 is a small increment of the variable k around the resonant 
value k , then we can expand the Hamiltonian, written in the new variables, in powers of 
K, and study the dynamics generated by the leading terms. We begin with the resonant 
part of the perturbation, 

SH^ = e 2 V}2\k,j) cos(4# + 25 + ip) + (terms with different n and m) 

— ^ £ 2 kgj cos(7 + ij; — ip Q ) + (terms with different n and m) . (3.11) 

Here, j is to be thought of as j(K = 0, J), with J arbitrary. Since this term in the 
perturbation depends only on 7 and not on 5, J is still a constant of the motion; for 
simplicity we fix its value at the resonant one, J . Then for the unperturbed Hamiltonian 
we obtain 

H (K, J ) = H (k ,j ) + \ O K 2 + 0(K 3 ) , (3.12) 



where 



_ d 2 H 2 d 2 H d 2 H 2 _ 8/j 2 fi 4 k 

n ~dk 2 ~ n ° + 2 djdk n ° m ° + ~w m ° ~ v/ i - fi 2 Jo + (1 - f, 2 Jo r/ 2 • (3 - 13) 
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So, to the lowest order in K and e, the K dynamics near the resonant torus is generated 
by 

tf loc ( 7 , K) = H (k , Jo ) + \ n K 2 + § e 2 k 2 j cos( 7 + V - V>o) , (3-14) 

which is the Hamiltonian of a well-known system, the non-linear pendulum. 

Let us now find the homoclinic orbits. If we call h Q = H (k ,j ), as before, and fix 
some value H loc = h loc for the local Hamiltonian, we can compute from Eq. (3.14) 

M 2 

Kc - h o - uf £2k o Jo cos (7 + i> - A) 



To simplify calculations, we set ip = i/j+tt. Then, for the homoclinic orbit, at the maximum 

VI 
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7 = n of the potential we must have K = 0, from which we find that h loc — h = ^§ e 2 kg j , 



and 



K = ±{fn- 1 e 2 k 2 o3o (1 + cos 7 )} 1/2 = ± n-hX 3o y /Z cos | . (3.15) 

Finally, we can compute the evolution of 7 by integrating the Hamilton equation 7 = 
dH loc /dK = VtK, and we obtain 



1/2 



7(£) = 4arctanexp j (^e 2 £lkl j j / £ j - tt . 



(3.16) 



We will use all of these results in the next section. 



IV. The Melnikov method 

We are now going to study the effect of the other terms of the perturbation 
SH'W ( 7 , 5; K, J) = e 2 J2 V«$ (kj) cos(n6 + m5 + ^ n J 

nm 

= e 2 £ V nm (*> J) cos 7 + S-^ + ^ nm ) , (4.1) 

nm ^ ' 

in particular of the terms with (n, m) 7^ (n , m ), on the dynamics of the integable system 
with Hamiltonian H loc , near the homoclinic orbit. We will use the Melnikov method to 
show that a stochastic layer forms in the vicinity of this destroyed separatrix, which acts 
as a seed for chaos. 

We thus want to show 13 that the Melnikov function fo the local K dynamics has a 
contribution 

M nm (^ nm ):=[ dt{H loc ,8H'M} (4.2) 
A.o. 

which is an oscillating function of ?/> nm and has therefore isolated, transverse zeroes, for 
some values of (n, m); here, the integral in (4.2) is taken over the homoclinic orbit, given 
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by in (3.16), 6(t) = Ujt + 0(e 2 ) from 6 = dH/dj', k = k Q + AK = k + 0(e), and 
j=j + O(e). 

In the calculation of the Melnikov integral we will keep only the leading order terms 
in e and K. Then, approximating the Poisson bracket 

2m — n 



{H loc , 6H'V>] ~ ji flK 2 , eVf cosQ 7 + ^ 5 - + Vw) } 

= 1 <g) (]fe , , ) sin (2 7 + * - + Vw) , (4-3) 

and using 7 = VtK, we get 

J — 7T 

= I ^ Km (*Wo) (^nm Sin V4m + 5 nm COS , 

where we have defined "4>' nm '■= ifJ nm — f V'o' ano - 

A n m - = J d7cos(^-7 + 5( 7 ) j 

„ Z" 7 " , . /n 2m — n . . . \ 
B nm : = J d 7 sin^-7 + 5( 7 )J, 

and 5(7) is obtained from (3.16) and S u^-t, 



(4.4) 



(4.5) 



5(7) = 77 In tan 



7 + 77 



-1/2 



(4.6) 



The function M nm in (4.4) is clearly oscillating whenever nV nm A nm and/or nV nT n B nm 
are not zero. We thus have to show that, for some (n, m) in the set {(2, 0), (2, ±2), (2, ±4), 
(4, 0), (4, -2), (4, ±4)}, A nm or B nm don't vanish. 

Perhaps surprisingly, the integrals in (4.5) can be calculated analytically, at least for 
half of the cases of interest, the ones with n = 2. To do this, notice that A nm and B nm 
are, respectively, the real and imaginary parts of the complex function 

i(2m-n)r?/2 



d7e 



tan 



7 + 77 



r n/2 

4e~ in7r/4 / dxe [nx (tanx) [{2m - n)ll/2 , 
Jo 



(4.7) 



where we have defined x := \ (7 + tt). Therefore, for n = 2, we have 17 



^2™ = -4i 



= -4i 



r n/2 -tt/2 

/ dx cos 2x tan 1 ^- 1 ^ x + i / dx sin 2x tan 1 ^" 1 ^ 



.7" 



i(m — 1)t]tv i(m — 1)7777 i(m — 1)7777 i(m — 1)7777 



— 2(m — 1)7777 



, (m — 1)7777 . (m — 1)7777 
sech h 1 cosech 



(4.8) 
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in other words, 

^2m = -2(m-l)?77rsech^ ^-U- , B 2m = -2(m - l)r/7r cosech ^ , (4.9) 

which don't vanish for any of the values of m listed above. 

V. Conclusions 

We have shown analytically that a RW universe, filled with a massive conformally coupled 
scalar field, has a chaotic behaviour, at least for a sufficiently small value of the mass, 
confirming the numerical results of Ref. 1. This is just the first step in our program to 
analyze the dynamics of RW models with various types of matter content, in order to relate 
this to numerical results, and make more precise statements about physical consequences. 
Among the purely theoretical aspects of this work, on the one hand, some possible gen- 
eralizations to cosmological models differing by the addition of a cosmological constant 
and/or other parameters are being currently studied, 18 as well as the possibility of adding 
a second matter field; on the other hand, we are working on extending the results obtained 
for the model in this paper to an estimation of the size of the stochastic region in phase 
space and of the time scale for the manifestation of the chaos we predict. The latter point 
is an important one since, as we mentioned in §2, this time scale must be smaller than the 
lifetime of the universe in our model. Although the simulations in Ref. 1 covered many 
cycles for the universe, it was stated there that chaos would show up on much smaller 
time scales, and therefore has, in principle, observable consequences; in Ref. 7, however, 
the authors argued qualitatively that they expect the time scale to be much larger than 
the lifetime of the universe, and used this to motivate the use of two different scalar fields 
in their model. We believe that the issue should be resolved by an actual, quantitative 
estimate of the time scale. 

We will conclude with some comments on possible avenues for future research, concern- 
ing other aspects of the relationship between chaos and the physics of the early universe. 
At a classical level already, chaos opens up ways of explaining the origin of the arrow of 
time, 19 and in our cosmological setting therefore, of improving our knowledge of cosmolog- 
ical statistical mechanics and thermodynamics. Furthermore, just as non-integrability and 
dynamical instability are, at the classical level, the causes of chaos, they are also, at the 
quantum level, the causes of the instability of quantum states and particles; one therefore 
expects, for example, classical chaos and semiclassical particle production to be related. 20 
At the semiclassical level, 21 one can prove the presence of decoherence and correlations, 
which cause the transition to the classical regime; stable and unstable scalar field quanta 
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are created during this process. The study of this passage, that leads from particle creation 
to chaos is, therefore, imperative. In a model very similar to ours, 22 a Lyapunov variable 
was used to define a growing entropy in the universe; this could be adapted, with minor 
changes, to the present model, and looks primising since, for example, growth of entropy, 
particle creation, decoherence and isotropization (a tendency to some kind of equilibrium) 
are related in semiclassical cosmo logical models. 23 

The relevance of chaos in the transition to the semiclassical limit of quantum cosmol- 
ogy, for RW models, was studied in Ref. 24. At the complete quantum level, let us recall 
only that RW models have been used to study the relation between the cosmological and 
thermodynamical arrows of time; 25 but many other results related to these subjects can be 
found in the literature. Ultimately, a very interesting goal for future work, and possibly a 
formidable task even in our simple RW model with chaotic behaviour, is to relate all these 
fundamental features of the universe in a unified theoretical framework. 
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